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ABSTRACT 

The spatial distribution of gas matter inside galaxy clusters is not completely smooth, 
but may host gas clumps associated with substructures. These overdense gas substructures are 
generally a source of unresolvedbias of X-ray observations towards high density gas, but their 
bright luminosity peaks may be resolved sources within the ICM, that deep X-ray exposures 
may be (already) capable to detect. In this paper we aim at investigating both features, using 
a set of high-resolution cosmological simulations with ENZO. First, we monitor how the bias 
by unresolved gas clumping may yield incorrect estimates of global cluster parameters and 
affects the measurements of baryon fractions by X-ray observations. We find that based on 
X-ray observations of narrow radial strips, it is difficult to recover the real baryon fraction to 
better than 10-20 percent uncertainty. Second, we investigated the possibility of observing 
bright X-ray clumps in the nearby Universe (z ^ 0.3). We produced simple mock X-ray ob- 
servations for several instruments (XMM, Suzaku and ROSAT) and extracted the statistics of 
potentially detectable bright clumps. Some of the brightest clumps predicted by simulations 
may already have been already detected in X- ray images with a large field of view. However, 
their small projected size makes it difficult to prove their existence based on X-ray morphol- 
ogy only. Preheating, AGN feedback and cosmic rays are found to have little impact on the 
statistical properties of gas clumps. 
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1 INTRODUCTION 

The properties of gas matter in more than two-thirds of the galaxy 
cluster volume are still largely unknown. In order to further im- 
prove our use of clusters as high-precision cosmological tools, 
it is therefore necessary t o gain insight into the t hermodynamics 
of their outer regions (see lEttori & MolendiluOllL and references 
therein). The surface brightness distributi on in clusters out skirts 
has been studied with ROSAT PSPC (e.g lEckert et all 1201 1 and 
references therein) and with Chandra (see lEttori & Balestral 2009. 
and references therein). A step forward has been the recent ob- 
servation of a handful of nearby clusters with the Japanese satel- 
lite Suzaku that, despite the relatively poor PSF and small field of 
view of its X-ray imaging spectrometer (XIS), benefits from the 
modest backgrou nd associated to its low-Ea rth orbit (see results 
on PKS0745-191. lGeorge et al.ll2009l: A2204. iReiprich et al.ll2009l : 
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A1795. iBautz et all 120091 : A14 13. iHoshino et all |201(A A1689, 
iKawaharada et alj20ld : Perseus JSimionescu et al. 201 1 ; Hydra A, 
ISato et al1l2012l : see also lAkamatsu & K awahara 2011). The first 
results from Suzaku indicate a flattening and sometimes even an 
inversion of the entropy profile moving outwards. The infall of 
cool gas from the large-scale structure might cause the assumption 
of hydrostatic equilibrium to break down in these regions, which 
could have important implications on cluster mass m easurements 
dRasia et al. 1120041 : lLau et aDl2009l : lBurns et alJ2010h . However, as 
a consequence of the small field of view (and of the large solid an- 
gle covered from the bright nearby clusters) observations along a 
few arbitrarily chosen directions often yield very different results. 



20101 : iMathews & Guol 1201 it iBattadia et all l2010l : IVazza et al. 
20l3;[Bode et al. 20121) . there are "simpler" physical reasons to ex- 



In addition to instrumental effects (e.g. Eckert et al. 2011), 
or no n -gravitational effects (e . g. iRoncarelli et alj|2006 \_ Lapi et al, 



pect that properties of the ICM derived close to ~ R200 may be 
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more complex than what is expected from idealized cluster mod- 
els. 

One effect is cluster triaxiality and asymmetry, which may 
cause variations in the ICM properties along directions, due to the 
presenc e of large-scale filaments in particular sectors of each clus- 
ter (e.g. IVazza et al J 1201 id : lEckert et al.ll2012h . or to the cluster- 
to-cluster varianc e related to the surrounding environment (e.g. 
iBurns et alJ2010h . 

A second mechanism related to simple gravitational physics is 
gas clumping, and its pos sible variations across different direction s 
from the cluster centre (Mathi esen et al.lll999l ; lNagai & Laull201lh . 

In general, gas clumping may constitute a source of uncer- 
tainty in the derivation of properties of galaxy clusters atmospheres 
since a significant part of detected photons may come from the 
most clumpy structures of the ICM, which may not be fully rep- 
resentative of the und e rlying large-scale d istribution of gas (e.g. 
iRoncarelli et alj 120061 : IVazza et all l201ld) . Indeed, a significant 
fraction of the gas mass in cluster outskirts may be i n the form of 
dense gas clumps, as suggested by recent simulations (iNagai & Laul 
201 H). In such clumps the emissivity of the gas is high, leading to an 
overestimation of the gas density if the assump tion of constant den- 
sity i n each shell i s mad e. The recent results of Nagai & Laul (1201 ll) 
and lEckert et all d2012h show that the treatment of gas clumping 
factor slightly improves the agreement between simulations and ob- 
served X-ray profiles. 

The gas clumping factor can also bias the derivation of the 
total hydr ostatic gas mass in g alaxy clusters at the ~ 10 per- 
cent level (Mat hiesen et a l. 1999), a nd it m ay as well bias the pro- 
jected temperature low jRasia et al] [2005b . Theoretical models of 
AGN feedback also suggest that overheated clumps of gas may 
lead to a more efficient d istribution of energy within cluster cores 
(iBirnboim & DekefcOllh . 

Given that the effect of gas clumping may be resolved or un- 
resolved by real X-ray telescopes depending on their effective res- 
olution and sensitivity, in this paper we aim at addressing both 
kind of effects, with an extensive analysis of a sample of massive 
(tsj 1O 15 M ) galaxy clusters simulated at high spatial and mass 
resolution (Sec|2]). First, we study in Sec B.ll the bias potentially 
present in cluster profiles derived without resolving the gas sub- 
structures. Second, in Sec l3.2l we derive the observable distribution 
of X-ray bright clumps, assuming realistic resolution and sensitiv- 
ity of several X-ray telescopes. We test the robustness of our results 
with changing resolution and with additional non-gravitational pro- 
cesses (radiative cooling, cosmic rays, AGN feedback) in Sec. 3. 3. 
Our discussion and conclusions are given in Sec|4] 



2 CLUSTER SIMULATIONS 

The simulations analysed in this work were produced with the 
Adaptive Mesh Refinement code ENZO 1.5, developed by the Lab- 
oratory for Computational Astrophysics at the University of Cali- 
fornia in San DiegoQ fe.g. lNorman et alJl2007l : ICollins et aDl2010l . 
and references therein). We simulated twenty galaxy clusters with 
masses in the range 6 • 10 14 ^ M/M < 3 • 10 15 , extracted from 
a total cosmic volume of Lbox ~ 480 Mpc/h. With the use of a 
nested grid approach we achieved high mass and spatial resolution 
in the region of cluster formation: radm = 6.76 • 1O 8 M for the 
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Figure 1. 2-dimensional power spectra of X-ray maps for three simulated 
clusters at z = (El, E3B and E15B, as in Fig0. Only one projection 
is considered. The long-dashed lines show the power spectra of the total 
X-ray image of each cluster, while the solid lines show the spectrum of 
each image after the average X-ray profile has been removed. Zero-padding 
has been considered to deal with the non-periodic domain of each image. 
The spectra are given in arbitrary code units, but the relative difference in 
normalization of each cluster spectrum is kept. The vertical grey line shows 
our filtering scale to extract X-ray clumps within each image. 



Table 1. Main characteristics of the simulated clusters at z = 0. Column 



1: identification number; 2: total virial mass (M v 



Mi 



DM ■ 



M gas ); 3: 



virial radius (R v ); 4:dynamical classification: RE=relaxing, ME=merging 
or MM=major merger; 5: approximate redshift of the last merger event. 



ID 




R v 


dyn. state 


redshif 




[1O 15 M ] 


[Mpc] 






El 


1.12 


2.67 


MM 


0.1 


E2 


1.12 


2.73 


ME 




E3A 


1.38 


2.82 


MM 


0.2 


E3B 


0.76 


2.31 


ME 




E4 


1.36 


2.80 


MM 


0.5 


E5A 


0.86 


2.39 


ME 
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0.66 


2.18 


ME 




E7 


0.65 


2.19 


ME 




Ell 


1.25 


2.72 


MM 


0.6 


E14 


1.00 


2.60 


RE 




E15A 


1.01 


2.63 


ME 




E15B 


0.80 


2.36 


RE 




E16A 


1.92 


3.14 


RE 




E16B 


1.90 


3.14 


MM 


0.2 


E18A 


1.91 


3.14 


MM 


0.8 


E18B 


1.37 


2.80 


MM 


0.5 


E18C 


0.60 


2.08 


MM 


0.3 


E21 


0.68 


2.18 


RE 




E26 


0.74 


2.27 


MM 


0.1 


E62 


1.00 


2.50 


MM 


0.9 



1 http://lca.ucsd.edu 



Properties of gas clumps and gas clumping factor in the ICM 3 



• 

relaxed * 


• 

post merger 


merging 


i 

© 

^> °* 

o 

relaxed * o 


* 

o « 



' o ^ * 

♦o ^ » 

^ 
o © 

post merger 


(i * 

merging V 





)0e-18 1.48e-15 7.41 e-1 5 3.11e-14 1.26e-13 5.00e-13 



Figure 2. Top panels: X-ray flux in the [0.5-2] keV (in [erg/(s • cm 2 )]) of three simulated clusters of our sample at z=0 (E15B-relax, El-post merger and 
E3B -merging). Bottom panels: X-ray flux of clumps identified by our procedure (also highlighted with white contours). The inner and outer projected area 
excluded from our analysis have been shadowed. The area shown within each panel is ~ 3 x 3 R200 for each object. 



DM particles and ~ 25 kpc/h in most of the cluster volume in- 
side the "AMR region" (i.e. y 2 — 3 P200 from the clu ster centre, 
see lVazza et al]l20ld : l\^^l201ld : lvazza et alJl20lial for further 
details). 

We assumed a concordance ACDM cosmology, with Qo = 
1.0, Q B = 0.0441, Qdm = 0.2139, 0\ = 0.742, Hubble parame- 
ter h = 0.72 and a normalization for the primordial density power 
spectrum of as = 0.8. Most of the runs we present in this work 
(Sec. 3. 1-3.2) neglect radiative cooling, star formation and AGN 
feedback processes. In Sec. 3. 3, however, we discuss additional runs 
where the following non-gravitational processes are modelled: ra- 
diative cooling, thermal feedback from AGN, and pressure feed- 
back from cosmic ray particles (CR) injected at cosmological shock 
waves. 

For consistency with our previous analysis on the same sam- 
ple of galaxy clusters (I Vazza et al]|201Cl . [2011al lcL we divided our 
sample in dynamical classes based on the total matter accretion 
history of each halo for z ^ 1.0. First, we monitored the time 
evolution of the DM+gas mass for every object inside the "AMR 
region" in the range 0.0 < z < 1.0. Considering a time lapse of 
At = 1 Gyr, "major merger" events are detected as total matter ac- 
cretion episode where M(£ + A£)/M(£) - 1 > 1/3. The systems 
with a lower accretion rate were further divided by measuring the 
ratio between the total kinetic energy of gas motions and the ther- 



mal energy inside the virial radius at z = 0, since this quantity pa- 
rame ter provides an indication of the dynam ical activity of a cluster 
(e.g. lTormen et "all ll997l :l Vazza et "all l2006h . Using this proxy, we 
defined as "merging" systems those objects that present an energy 
ratio > 0.4, but did not experienced a major merger in their past 
(e.g. they show evidence of ongoing accretion with a companion 
of comparable size, but the cores of the two systems did not en- 
counter yet). The remaining systems were classified as "relaxed". 
According to the above classification scheme, our sample presents 
4 relaxed objects, 6 merging objects and 10 post-merger objects. 

Based on our further analysis of this sample, this classifica- 
tion actually mirrors a different level of dynamical activity in the 
subgroups, i.e. rel axed systems on average host weaker shocks 
(I Vazza et aD l201(t) . they are characterized b y a lowest turbulent 
to thermal energy ratio dVazza et al.1 1201 lal) , and they are char- 
acterized b y the smallest amount of azimuthal s catte r in the gas 
propert ies (IVazza et al]l201ld : lickert et al.ll2012h . In I Vazza et aD 
(2011c) the same sample was also divided based on the analysis of 
the power ratios from the multi-pole decomposition of the X-ray 
surface brig htness images (P3/P0) , and the centroid shift (w), as 
described bv lBohringer et all (2010). These morphological param- 
eters of projected X-ray emission maps were measured inside the 
innermost projected 1 Mpc 2 . This leads to decompose our sam- 
ple into 9 "non-cool-core-like" (NCC) systems, and 1 1 "cool-core- 
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like" systems ( CC), once that fiduc ial thresholds for the two pa- 
rameters (as in ICassano et al are set. We report that (with 
only one exception) the NCC-like class here al most perfect l y over - 
lap with the class of post-merger systems of I Vazza et alJ (120101) . 
while the CC-like class contains the relaxed and merging classes of 
IVazza etalJbOlOh . 

Table 1 lists of all simulated clusters, along with their main 
parameters measured at z = 0. All objects of the sample have a 
final total mass > 6 • 1O 14 M , 12 of them having a total mass 
> 1O 15 M . In the last column, we give the classification of the 
dynamical state of each cluster at z = 0, and the estimated epoch 
of the last major merger event (for post-merger systems). 

2.1 X-ray emission 

We simulated the X-ray flux (Sx) from our clusters, assuming a 
single temperature plasma in ionization equili brium within each 3D 
cell. We use the APEC emission model (e.g. ISmith et al .11200 lb to 
compute the cooling function A(T, Z) (where T is the temperature 
and Z the gas metallicity) inside a given energy band, including 
continuum and line emission. For each cell in the simulation we 
assume a constant metallicity of Z/Zq = 0.2 (which is a good ap- 
proximation of the observ ed metal abundance in cluster outskirts, 
iLeccardi & Mol endi 2008). While line cooling may be to first ap- 
proximation not very relevant for the global description of the hot 
ICM phase (T ~ 10 8 K), it may become significant for the emis- 
sion from clumps, because their virial temperature can be lower 
than that of the host cluster by a factor ~ 10. Once the metallicity 
and the energy band are specified, we compute for each cell the X- 
ray luminosity, Sx = nHn e A(T, Z)dV, where uh and n e are the 
number density of hydrogen and electrons, respectively, and dV is 
the volume of the cell. 



2.2 Definition of gas clumps and gas clumping factor 

Although the notion of clumps and of the gas clumping factor is 
often used in the recent literature, an unambiguous definition of 
this is non-trivial 0. In this work we distinguish between resolved 
gas clumps (detected with a filtering of the simulated X-ray maps) 
and unresolved gas clumping, which we consider as an unavoid- 
able source of bias in the derivation of global cluster parameters 
from radial profiles. On the theoretical point of view, gas clumps 
represent the peaks in the distribution of the gas clumping factor, 
and can be identified as single "blobs" seen in projection on the 
cluster atmosphere if they are bright-enough to be detected. While 
resolved gas clumps are detected with a 2-dimensional filtering 
of the simulated X-ray images (based on their brightness contrast 
with the smooth X-ray cluster emission), the gas clumping factor 
is usually estimated in the literature within radial shells from the 
cluster centre. The two approaches are not fully equivalent, and in 
this paper we address both, showing that in simple non-radiative 
runs they are closely related phenomena, and present similar 
dependence on the cluster dynamical state. 



2 While this article was under review, Zhuravleva et al.l {2012) published a 
work in which they also characterize the level of inhomongeities in the sim- 
ulated ICM, based on the radial median value of gas density and pressure. 
They study the impact of cluster triaxiality and gas clumping in the deriva- 
tion of global cluster properites such like the gas mass and the gas mass 
fraction, reporting conclusions in qualitative agreement with our work. 



2. 2. 1 Resolved gas clumps 

We identify bright gas clumps in our cluster sample by post- 
processing 2-dimensional mock observations of our clusters. We do 
not consider instrumental effects of real observations (e.g. degrad- 
ing spatial resolution at the edge of the field of view of observation), 
in order to provide the estimate of the theoretical maximum amount 
of bright gas clumps in simulations. Instrumental effects depends 
on the specific features of the different telescope, and are expected 
to further reduce the rate of detection of such X-ray clumps in real 
observations. In this section we show our simple technique to pre- 
liminary extract gas substructures in our projected maps, using a 
rather large scale (300 kpc/h). "X-ray bright gas clumps", in our 
terminology, correspond to the observable part of these substruc- 
tures, namely their small « 50 kpc/h) compact core, which may 
be detected within the host cluster atmosphere according the effec- 
tive resolution of observations (Sec l3.2l ). 

Based on the literature (e.g. iDolag et al.l2005l . l2009h . the most 
massive substructures in the ICM have a (3 -dimensional) linear 
scale smaller than < 300 — 500 kpc/h. We investigated the typical 
projected size of gas substructures by computing the 2-dimensional 
power spectra of Sx for our cluster images. In order to remove the 
signal from the large-scale gas atmosphere we subtracted the aver- 
age 2-dimensional cluster profile from each map. We also applied 
a zero-pa dding to take into ac count the non-periodicity of the do- 
main (see 1 Vazza et al.ll201 lal and references therein). The results 
are shown in FigQ] for three representative clusters of the sample 
at z = 0: the relaxed cluster E15B, the post-merger cluster El 
and the merging system E3B -merging (see also Fig0. The long- 
dashed profiles show the power spectra of the X-ray image of each 
cluster, while the solid lines show the spectra after the average 2- 
dimensional profile of each image has been removed. The power 
spectra show that most of the residual substructures in X-ray are 
characterized by a spatial frequency of k > 10 — 20, correspond- 
ing to typical spatial sca les lp < 300 — 6 00 kpc/h, similar to 
three-dimensional results (IDolag et al.l2QQ9h . A dependence on the 
dynamical state of the host cluster is also visible from the power 
spectra: the post-merger and the merging clusters have residual X- 
ray emission with more power also on smaller scales, suggesting 
the presence of enhanced small-scale structures in such systems. 
We will investigate this issue in more detail in the next sections. 

To study the fluctuations of the X-ray flux as a result of the 
gas clumps we compute maps of residuals with respect to the 
X-ray emission smoothed over the scale lo. The map of clumps 
is then computed with all pixels from the map of the residuals, 
where the condition Sx/Sx, smooth > a is satisfied. By impos- 
ing lo = 300 kpc/h (~ 12 cells) and a = 2, all evident gas 
substructures in the projected X-ray images are captured by the al- 
gorithm. We verified that the adoption of a larger or a smaller value 
of lo by a factor ~ 2 does not affect our final statistics (Sec l3.2l ) 
in a significant way. In Figure [2 we show the projected X-ray flux 
in the [0.5-2] keV energy range from three representative clusters 
of the sample (El -post merger, E3B -merging and E15B -relaxed, in 
the top panels) at z = 0, and the corresponding maps of clumps 
detected by our filtering procedure (lower panels). The visual in- 
spection of the maps show that our filtering procedure efficiently 
removes large-scale filaments around each cluster, and identifies 
blob-like features in the projected X-ray map. The relaxed system 
shows evidence of enhanced gas clumping along its major axis, 
which is aligned with the surrounding large-scale filaments. In- 
deed, although this system is a relaxed one based on its X-ray mor- 
phology within R200 1 2 (based on the measure of X-ray morpho- 
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logical parameters, Cassano et al. 2010, and to its very low value 
of gas turbulence, |Vazza et al. 2011 J), its large-scale environment 
shows ongoing filamentary accretions and small satellites, which is 
a rather common feature even for our relaxed cluster outside R500 . 
The presence of clumps in the other two systems is more evident, 
and they have a more symmetric distributions. We will quantify the 
differences of the gas clumping factor and of the distribution of 
bright clumps in these systems in the following sections. 

2.2.2 The gas clumping factor 

The definition of the gas clumping factor, C p , follows from the 
computation of the cluster density profile averaged within radial 
shell: 



all clusters z = 



10 



C P (R) ^ 



1 J n P 2 (R)dn 



(1) 



where at each radial bin from the cluster centre (which we define 
according to the peak of the gas mass density), R, we compute 
the angular average within radial shells (with constant width of 
1 cell at the highest AMR level, equivalent to 25 kpc/h) from 
the centre of clusters out to « 1.5 R200. This definition of the 
gas clumping factor is often used to interpret observed depar- 
tures from the smooth gas density, suggested by some observations 
dSimionescu et alJl2QllL iNagai & Laull201lb . However, averaging 
within spherical shells to compute C P (R) is a procedure prone to 
errors in the cases of mergers or asymmetries in the cluster atmo- 
sphere, because these phenomena break the spherical symmetry in 
the cluster and the a spherical average might not be a good ap- 
proximation. A-priori, the presence of a large gas clumping fac- 
tor (measured as in Eq[T} and the increased presence of dense gas 
clumps may be not associated phenomena. However, as we will see 
in the following Sections, a larger statistics of gas clumps and the 
presence of large-scale asymmetries are closely related phenomena, 
which are regularly met in perturbed galaxy clusters. Indeed, the re- 
moval of the radial average of the X-ray profile or the filtering of 
^ 300 kpc/h structures in X-ray images in FigOvisibly highlights 
the same structures, because all bright gas clumps are found within 
the sectors where the largest departure from the azimuthal profile 
are present (i.e. due to the presence of large-scale filaments). In the 
next Sections, we will explore the trend with the cluster dynami- 
cal state (and other cluster parameters) of the gas clumping factor 
and of the distribution of bright clumps, and we will suggest likely 
observational implications of both complementary by-products of 
substructures within clusters. 



3 RESULTS 

3.1 Gas clumping factor from large-scale structures 

The average clumping factor of gas in clusters is a parameter 
adopted in the theoretical interpretation of d e-projected profiles 
of gas mass, tempera t ure and entropy (e.g. lUrban et al] 1201 ll : 



ISimionescu et al.ll201 ll : lEckert et al.ll2012l : I Walker -et al.ll2012h 

For each simulated cluster, we compute the profile of the gas 
clumping factor following Sec l2.2.2l 

Figure [3] shows the average profile of C P (R) for the three dy- 
namical classes into which our sample can be divided. The trends 
for all classes are qualitatively similar, with the average clumping 
factor C p < 2 across most of the cluster volume, and increasing to 
larger values (C p ^3 — 5) outside R200 . Post-merger and merging 
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Figure 3. Profile of azimuthally averaged gas clumping factors for the sam- 
ple of simulated clusters at z = (the error bars show the la deviation). 
The different colours represent the three dynamical classes in which we 
divide our sample (four relaxed, six merging and ten post-merger clusters). 



systems present a systematically larger average value of gas clump- 
ing factor at all radii, and also a larger variance with in the subset. 
These findings are consistent with the recent results of lNagai & Laul 
(1201 11) . By extracting the distribution of bright gas clumps in pro- 
jected cluster images, we will see in Sec l3.2l that both the radial 
distribution of clumps, and the trend of their number with the clus- 
ter dynamical state present identical behaviour of the gas clumping 
factor measured here. This again support the idea that the two phe- 
nomena are closely associated mechanism, following the injection 
of large-scale substructures in the ICM. 

Next, we calculate the average profiles of gas clumping factor 
for different bins of gas over-density and temperature (limiting to 
the relaxed and the merging subset of clusters). This way we can 
determine if the gas clumping factor affects all phases of the ICM 
equally. 

In Fig|4]we show the average profiles of gas clumping factor 
calculated within different bins of gas over-density (£p C r,b < 50, 
50 < 5p cr , h < 10 2 , 10 2 ^ 5p cr , h < 10 3 and 5p cr , h ^ 10 3 , 
where 5p cv ^ is p/(/bPcr), with fh cosmic baryon fraction and p cr 
the cosmological critical density) and gas temperature (T < 10 6 
K, 10 6 < T < 10 7 K, 10 7 ^ T < 10 8 K and T ^ 10 8 K) at 
z = 0. 

We want to have a characterization of the environment asso- 
ciated with the most significant source of gas clumping (e.g. gas 
substructures) in a way which is unbiased by their presence. This 
is non-trivial, because if the local overdensity is computed within a 
scale smaller than the typical size of substructures « 300 kpc/h), 
the gas density of a substructure will bias the estimate of local over- 
density high. For this reason, the local gas overdensity and gas 
temperature in Fig|4]have bee computed for a much larger scale, 
1 Mpc/hB In FigH we also show the gas mass fraction of each 



We checked that the use of the DM over density statistics, in this case, 
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Figure 4. Average profiles of gas clumping factor and gas mass distribution for different phases across the whole cluster sample at z = 0, for relaxed clusters 
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phase inside the radial shells. The gas phases show evidence that 
the increased clumping factor of merging clusters is related to the 
increased clumping factor of the low-density phase, at all radii, 
even if the mass of this component is not dominant in merging sys- 
tems. An increased clumping factor is also statistically found in the 
cold « 10 6 K) and intermediate 10 6 < T < 10 7 phase of gas 
in merging systems. We remark that although the increase in this 
phase is associated to an enhanced presence of X-ray detectable 
clumps (Sec l3.2h it cannot be directly detected in X-ray, due to its 
low emission temperature. 

The gas phase characterized by 5p CY: h < 50 and T ^ 

10 6 - 10 7 is typical of large-scale fil aments toolag etal.ll2008l : 
IVazza et al.ll2009l : llapichino et al.ll201ll) . While filaments only host 
a small share of the gas mass inside the virial radius, they can dom- 
inate the clumping factor within the entire volume of the merging 
systems. This is because their gas matter content contains signifi- 
cantly more clumpy matter compared to the ICM onto which they 
accrete. This suggests that a significant amount of gas substructure 
(e.g. including massive and bright X-ray clumps, as in Sec j3.21 ) can 
originate from low-density regions, and not necessarily from the 
hot and dense phase of the ICM. This is reasonable in the frame- 
work of hierarchical structure formation, since matter inside fila- 
ments had almost the enti re cosmic evolution to develop substruc- 
tures jEinasto et al]|201lh . Contrary to most of the matter inside 
galaxy clusters, large-scale filaments have never undergone a phase 
of violent relaxation and efficient mixing. Therefore, even at late 
cosmic epochs they can supply very clumpy material. This sug- 
gests that clusters in an early merging phase are characterized by 
the largest amount of gas clumping factor at all radii, and that this 
clumpy materials is generally associated with substructures coming 
from filaments and from the cluster outskirts. Filaments between 
massive galaxy cl usters are indeed fo und to anticipate major merg- 
ers in simulations dVazza et al.ll2011ah . starting the injection of tur- 
bulence in the ICM already some ~ Gyr before the cores collision. 



We now investigate how gas clumping can affect X-ray ob- 
servations. For this purpose, we compare in Fig|5] the real pro- 
files of clumping factor, DM, gas density and baryon fraction for 
the same clusters as in FigQ with profiles derived from the en- 
tire cluster volume or within smaller volumes (i.e. thin slices along 
the line of sight) along perpendicular sectors of the same systems. 
All profiles are drawn from the centre of total (gas+DM) mass of 
each system. These selected volumes are meant to broadly com- 
pare with the selection of sectors chosen in recent deep X-ray ob- 
servations of clusters ( Simiones cu et al ] |201ll : lUrban etalJlioTTI : 
iHumphrev et al.ll2012h . Therefore, in this case we studied the ra- 
dial profiles insi de rectangular sectors of width ~ 500 kpc , similar 
to observations ( Simionescu et al .l201ll : IUrbanetal.l201lh . and the 
largest available line of sight (~ 8 Mpc). The projected total vol- 
ume sampled by these sectors is ~ 20 percent of the total cluster 
volume, and it is interesting to study how representative the infor- 
mation is that can be derived from them. 

For comparison with the observationally derived baryon 
fraction, we comp ute the best-fit model for NFW profiles 
dNavarro et al.| [l996) using the CURVEFIT task in IDL for the total 
volume, and for each sector independently. The radial range used to 
compute the best-fit to the NFW profile is 0.02 R 2 oo < R < R200. 



yield very similar results to those obtained using the gas over density. How- 
ever, we focus on the latter in this work, since this quantity can be more 
easily related to observations 
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In the relaxed system (first panels) there is little gas clump- 
ing out to a radius 0.8i?2oo- The average profile of DM is overall 
well modelled by the NFW profile everywhere up to ~ R200 . The 
profiles inside smaller sectors are also well described by a NFW 
profile, yet for > 0.6i?2oo the NFW profiles within each sector 
start to deviate from one another. At large radii the best-fit profiles 
are found to be both under- or over-estimating the real DM profile, 
due to large-scale asymmetries in the cluster atmosphere (e.g. left 
panel of Fig0. Asymmetries in this relaxed cluster also cause a 
si milar azimuthal scat ter of gas profiles, similar to what was found 
in lVazza et all (12011 ch . 

The departures from the NFW profile and the differences be- 
tween sectors are significantly larger in the post-merger and the 
merging systems (second and third columns of Fig[5]). We find that 
the best-fit NFW profiles can differ significantly between sectors, 
with a relative difference of up to ~ 30 — 40 percent at R200 in the 
merging system caused by the infalling second cluster (see right 
panel of Fig0. In addition to azimuthal variations, also large dif- 
ferences (> 40 percent) between the true DM profile and the best- 
fit NFW profile are found along some sectors. 

These uncertainties in the "true" gas/DM mass within se- 
lected sectors bias the estimated enclosed baryon fraction. This 
is an important effect which may provide the solution to un- 
derstand the variety of recent results provided by observations 
( Simionescu et al ] |201ll : lUrban et al.ll20lTl : iHumphrev et aDl2012l : 
IWalker et al.ll2012h . 

In Fig(5j we present the radial trend of the enclosed baryon 
fraction, Yg as (< R) = /bar(< R)/ fb for each cluster and sector 
(where fh is the cosmic baryon fraction), following three different 
approaches: 

• we measure the "true" baryon fraction inside spherical shells 
(or portion of spherical shells, for the narrow sectors); 

• we measure the "X-ray-like" baryon fraction given by 
(/ 47rR 2 p(R) 2 dR) - 5 / J 47rR 2 p(R) NFW dR, where p N Fw is the 
profile given by the best-fit NFW profile for DM (for the total vol- 
ume or for each sector separately); 

• we measure the "clumpy" baryon fraction, 
(/ 47vR 2 p(R) 2 dR/ J 47iR 2 p(R)l M dR) - 5 , where P (R)um 
is the true DM density). 

This "clumpy" estimate of Y gas obviously has no correspond- 
ing observable, since it involves the clumping factor of the DM dis- 
tribution. However, this measure is useful because it is more con- 
nected to the intrinsic baryon fraction of the massive substructures 
within galaxy clusters. 

The radial trend of the true baryon fractio n for the total vol- 
ume is in line with t he result of obser vations (I Allen et all 2002; 
IVikhlinin et"ai]l2006l : lEttori et al.ll2009h . with Y s increasing from 
very low central values up to the cosmic value close to R200 . The 
total enclosed baryon fraction inside < 1.2 R200 is ±5 percent of 
the cosmic value, yet larger departures (±10 — 20 percent) can be 
found inside the narrow sectors. 

The "X-ray-like" measurement of the enclosed baryon frac- 
tion is subject to a larger azimuthal scatter, and shows larger differ- 
ences also with respect to the true enclosed baryon fraction. In our 
data- sample, it rarely provides an agreement better than ±10 per- 
cent with the true enclosed baryon fraction (considering the whole 
cluster or the sectors). Based on our sample, we estimate that on 
average in ^ 25 percent of cases the baryon fraction measured 
within sectors by assuming an underlying NFW profile overesti- 
mates the true baryon fraction by a ~ 10 percent, while in ~ 50 



percent of cases it underestimates the cosmic baryon fraction by a 
slightly larger amount (~ 10 — 20 percent). In the remaining cases 
the agreement at R200 is of the order of ^ 5 percent. 

In all cases the main factor for significant differences between 
the X-ray-derived baryon fraction and the true one is the triaxi- 
ality of the DM matter distribution within the cluster rather than 
enhanced gas clumping. As an example, in the relaxed cluster of 
Fig(5]the "X-ray" baryon fraction is underestimated in 2 sectors, it 
is overestimated in 1 sector and it almost follows the true one in the 
last sector. This trend is mainly driven by the difference between 
the true DM profiles and the best-fit NFW profiles within sectors 
for > 0.4 i?2oo, rather than dramatic differences in the gas clump- 
ing factor. 

In general, a significant gas clumping factor (C p ~ 2 — 3) in 
the cluster outskirts can still lead to fairly accurate baryon fractions 
provided that this clumping is associated with DM substructures or 
DM filaments. In a realistic observation, the presence of a large- 
scale filament towards the cluster centre would still yield a reason- 
able best-fit NFW profile, even if not necessarily in agreement with 
the global NFW profile of the cluster. As a result, the measurement 
of a baryon fraction close to the cosmic value along one sector does 
not imply the absence of gas clumping. On the other hand, the pres- 
ence of significant large-scale asymmetries can be responsible for 
strong differences between the baryon fraction measured in X-rays 
and the true baryon fraction, even without substantial gas clumping 
factor. To summarize, the result of our simulations is that in many 
realistic configurations obtained in large-scale structures, the gas 
clumping factor inside clusters and the enclosed baryon fraction 
may provide even unrelated information. 

Also the radial range selected to fit the NFW profile to the 
observed data can significantly affect the estimate of the enclosed 
baryon fraction. Given the radial gas and DM patterns found in our 
simulated clusters, all radii outside > 0.4i?2oo may be character- 
ized by equally positive or negative fluctuations compared to the 
NFW profile. Adopting a different outer radius to fit a NFW profile 
to the available observations may therefore cause a small but not 
negligible (~ 5 — 10 percent in the normalization at the outer radii) 
source of uncertainty in the final estimate of the enclosed baryon 
fraction. 

These findings might explain why the enclosed baryon frac- 
tion is estimated to be larger tha n the cosmic baryon f raction at 
~ R200, for some nearby clusters dSimionescu et al.l201ll). smaller 
than the cosmic value in some other cases jKawaharad a^etaD 
120101) . and compatible with the cosmic value in th e remaining 
cases (e.g. Humphrey et al ] boil I Walke r et al ]|20l3)- These ob- 
servations could generally only use a few selected narrow sectors, 
similar to the procedure we mimic here, and the chance that the 
profiles derived in this way are prone to a ±10 — 20 percent uncer- 
tainty in the baryon fraction appears high. 

We conclude by noting that the "clumpy" baryon fraction in 
all systems is systematically lower by ^ 30 — 40 percent com- 
pared to the baryon fraction of the whole cluster. This suggest that 
the clumpiest regions of the ICM are baryon-poor compared to the 
average ICM. The most likely reaso n for that is that clumps are usu- 
ally a ssociated with substructures dTormen et a l. 2003; Dolag ^taD 
2009), that have been subjected to ram-pressure stripping. Most 
of these satellites should indeed lose a sizable part of their gas 
atmosphere while orbiting within the main halo, ending up with 
a baryon fraction lower than the cosmic value. Recent observa- 
tions of gas-depleted galaxies identified in 300 nearby groups with 
the Sloan Digital S ky Survey also seem to support this scenario 
(Zhang et al. 2012). This mechanism is suggested by numerical 
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Figure 8. Different distribution functions of clumps detected in our simu- 
lated X-ray maps at z=0.1 and assuming >Sx,iow — 2 • 10 -15 erg/(s • cm 2 ) 
and an effective spatial resolution of w 10" (to mimic a deep survey with 
XMM). The continuous lines are for the differential distribution, the dashed 
lines are for the cumulative distributions (only the total is shown). The top 
panel show the different luminosity functions for three different bin in the 
projected radii, the bottom panel shows the decompositions of the dataset 
into three total mass range. The dot-dash lines show the ratio of objects 
found in each class, compared to the average population. 



3.2 Properties of gas clumps 

From our simulations we can compute the number of resolved gas 
clumps in simulated X-ray images. Based on the clump detection 
scheme outlined in Sec l2.2.U we produced catalogues of X-ray 
bright clumps in our images, and we compared the statistics of dif- 
ferent detection thresholds and spatial resolution, broadly mimick- 
ing the realistic performances of current X-ray satellites. We con- 
sidered a range of redshifts (z — C0, z — 0.1 and z — 0.3), and 
accounted for the effect of cosmological dimming (oc (1 + z) 4 ) 
and of the reduced linear resolution. We did not generate photons 
event files for each specific observational setup, but we rather as- 
sumed a reference sensitivity, 5x,iow, for each instrument, which 
was derived from the best cases in the literature. The production 
of more realistic mock X-ray observations from simulated images, 
calibrated for each instruments, inv olves a number of technicali- 
ties and problems in m odelling (e.g. lRasia et al.ll2005l ; iHeinz et al.l 
l20ld : lBiffietai1l2012h . which are well beyond the goal of this pa- 
per. Our aim, instead, is to assess the maximum amount of bright 
clumps that can be imaged within the most ideal observational con- 
dition for each instrument. More realistic observational conditions 
(e.g. a decreasing sensitivity of real instruments moving outwards 
of the field of view, etc) are expected to further reduce the real 
detection ratio of such clumps (if any), and therefore the results 
presented in this Section should be considered as upper-limits. 

Our set of synthetic observations consists of 60 images (three 
projections per cluster) for each redshift and instrument. In Fig- 
ure [6] we present the luminosity function of X-ray clumps de- 
tected in our maps, assuming three different instrumental setups 
(we consider here the maximum luminosity per pixel for each 
clump), representative of realistic X-ray exposures: 10'Vpixel and 
Sx,iow = 2 • 10 _15 erg/(s • cm 2 ) for XMM-like observations, 
30"/pixel and Sx,iow = 3 • 10- 14 erg/(s • cm 2 ) for ROSAT-hke 
observations and 120"/pixel and Sx,iow = 10 _13 erg/(s • cm 2 ) 
for Suzaku-like exposures (for simplicity we considered the [0.5-2] 
keV energy range in all cases). In the case of the XMM "observa- 
tion" at the two lowest redshifts, it should be noted that the avail- 
able spatial resolution is significantly better than the one probed 
by our simulated images. In this case, a simple re-binning of our 
original data to a higher resolution (without the addition of higher 
resolution modes) has been adopted. 

Based on our results, a significant number of clumps can be 
detected in deep exposures with the three instruments: ^20 — 30 
per cluster with ROSAT/Suzaku, and ~ 10 2 per cluster with XMM. 
These estimates are based on observations that sample the whole 
cluster atmosphere (1.2 x 1.2 ifeoo). The evolution with redshift 
reduces the number of detectable clumps by ~ 2 — 3, mainly due 
to the effect of the degrading spatial resolution with increasing dis- 
tance (and not to significant cosmological evolution). At z = 0.3, 
only XMM observations would still detect some clumps (~ 30 per 
cluster). The largest observable flux from clumps, based on these 
results, is ~ 5 • 10 _12 erg/(s • cm 2 ) for Suzaku and ROSAT and 
- 10" 12 erg/(s • cm 2 ) for XMM0. 

In the following we restrict ourselves to the fiducial case of 



simulations (iTormen et al.ll2003l : iDolag et all 2009h, and it is very 
robust against numerical details ( e.g. [Y azza et al. 2011b) or the 
adopted physical implementations (IDolag et al.ll2009h . However, it 
seems very hard to probe observationally this result, because these 
DM sub-halos need high resolution data to resolve them both in 
their gas component with X-ray and in the total mass distribution 
through, e.g., gravitational lensing techniques. 



4 More exactly, we adopted z = 0.023, corresponding to the luminosity 
distance of the COMA cluster, w 100 Mpc. 

5 We note that the maximum observable flux here is not completely an 
intrinsic property of the clumps, but it also slightly depends on the post- 
processing extraction of clumps in the projected maps at different resolu- 
tion, as in Sec |2.2.1| For large PSF, indeed, the blending of close clumps 
can produce a boost in the detectable X-ray flux within the PSF. 
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Figure 7. Differential distribution functions of clumps detected in our simulated X-ray maps at z=0.1 and assuming *Sx,iow = 2 • 10 ~~ 15 erg/ (s • cm 2 ) and an 
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population. 



z — 0.1 observations with XMM, and we study the statistical dis- 
tributions of clumps in more detail. 

In Figure |7] we show differential distribution functions for the 
various parameters available for each detected clump: a) maximum 
X-ray flux of each identified clump; b) total X-ray flux on the 
pixel scale; c) projected clump radius (assuming a full-width-half- 
maximum of the distribution of SxEb; d) projected distance from 
the cluster centre; e) average projected temperature; f) temperature 
contrast around the clump (estimated by dividing the temperature 
of the clump by the average surrounding temperature of the ICM, 
calculated as the average after the exclusion of the clump tempera- 
ture). 

The luminosity distributions of clumps confirm the result of 



6 The projected clump radius is measured here with an extrapolation, by as- 
suming that the observed flux from clumps originates from a Gaussian dis- 
tribution, for which we compute the full-width-half-maximum. This choice 
is motivated by the fact that in most cases our detected clumps are very 
close to the resolution limit of our mock observations, and only with some 
extrapolation we can avoid a too coarse binning of this data. For this rea- 
son, our measure of the clump radius must be considered only a first order 
estimate, for which higher resolution is required in order to achieve a better 
accuracy. 



the previous section, and show that post-merger and merging clus- 
ters are characterized by a larger number of detectable clumps. 
They also host on average clumps with a ~ 2 — 3 times larger X-ray 
flux. Perturbed clusters have a factor ^2 — 4 excess of detectable 
clumps at all radii compared to relaxed clusters, and they host larger 
clumps. Clumps in post-merger systems are on average by factor 
~ 2 smaller than in merging systems, likely as an effect of strip- 
ping and disruption of clumps during the main merger phase. The 
distributions of temperatures and temperature contrasts show that 
the projected clumps are always colder than the surrounding ICM. 
This is because the innermost dense region of surviving clumps is 
shielded from the surrounding ICM, and retains its lower virial tem- 
perature for a few dynamical times dTormen et al .1 12003). For most 
of our detected clumps, the estimated FWHM is ^ 50 kpc/h. This 
is close to the cell size in our runs, and raises the problem that the 
innermost regions of our clumps are probably not yet converged 
with respect to resolution. We will revisit this issue again in the 
next section. 

In Sec(8] we investigate the dependence of the luminosity 
function of clumps on the projected radial distance from the centre 
(top panel) and on the total gas mass. We found that for all luminos- 
ity bins the intermediate annulus (0.6 ^ R/R200 < 1.2) contains 
> 70 percent of detectable clumps, even if its projected volume is 
about the same as that of the most external annulus. On the other 
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hand we detect no significant evolution with the host cluster mass. 
This is reasonable since we neglect radiative cooling or other pro- 
cesses than can break the self- similarity of our clusters. 

Given the relatively large number of detectable clumps pre- 
dicted by our simulations, one may wonder whether existing ob- 
servations with XMM, Suzaku, ROSAT or Chandra can already be 
compared with our results. We will come back to this point in Sec|4] 

3.3 Tests with additional physics and at higher resolution 

Several physical and numerical effects may affect the amount 
of gas clumping factor measured in simulations. For instance, 



Nagai &Laul(l201 ]]) recently showed that the number of clumps in 
simulations with radiative cooling and feedback from supernovae 
and star formation is significantly larger than in non-radiative runs. 
However, clumps can cool so efficiently that they can even tually 
cool below X-ray emitting temperatures (Na gai & Laull201 lb . Ra- 
diative cooling leads to the formation of more concentrated high- 
density clumps in the ICM, while feedback from SN and AGN may 
tend to wash them out by preventing the cooling catastrophe and 
providing the gas of clumps more thermal energy. The result of this 
competition between radiative losses and energy input via feedback 
may depend on the details of implementation. Also, the presence 
of cosmic rays accelerated at cosmological shocks may yield a dif- 
ferent compressibility of the ICM in the cluster outsk irts, which 
can c hange significantly the amount of clumping there (I Vazza et all 
l2012h . 

In this section, we assess the uncertainty in our results, using 
a set of cluster re- simulations with different setups. Cluster runs 
such as the one we discussed in the previous sections are fairly 
expensive in terms of CPU time (~ 3 — 4 • 10 4 CPU hours for 
each cluster), given the number of high-resolution cells generated 
during the simulation. In order to monitor the effects of different 
se tups we thus opted for a set of smaller clusters, already explored 
m Vazza (2011b). We chose a pair of galaxy clusters with final mass 
of ~ 2 — 3 • 1O 14 M and two very different dynamical states (one 
post-merger and one relaxed cluster), and studied in detail the radial 
distribution of gas clumping factor in each of them. These clusters 
were re- simulated with the following numerical setups: a) a sim- 
ple non-radiative run, with maximum resolution of 25 kpc/h; b) a 
run with pure radiative cooling, same peak resolution; c) a run with 
thermal feedback from a single, central AGN starting at z = 1, with 
a power of Wj e t ~ 10 44 erg/s for each injection, with same peak 
resolution; d) a run with a uniform pre-heating of 100 keVcm 2 (re- 
leased at z = 10) followed by a phase of thermal feedback from a 
central AGN starting at z = 1 and a power of Wj e t ~ 10 43 erg/s 
for each injection ; e) a run with cosmic-ray injection at shocks, re- 
duced thermalization and pressure feedback, same peak resolution; 
f) a simple non-radiative run at a higher res olution , 12 kpc /h. 

The cases c) and d) are taken from Vazza J2011bh . where 
we modelled the effect s of a uniform pre -heating background of 
cosmic gas (following lYounger & Brvanl l2007h and of intermit- 
tent thermal AGN activity, modelled as bipolar outputs of over- 
pressurized gas around the peak of cluster gas density. These runs 
were designe d to fit the average therm al properties of nearby galaxy 
clusters (e.g. ICavagnolo et al] l2009h and to quench catastrophic 
cooling in radiative cluster simulations. We note that the amount 
of pre-heating and thermal energy release from the "AGN region" 
were imposed by construction (in order to allow a parameter space 
exploration) and are not self-consistently derived from a run-time 
model of super-massive black-holes evolution. In this implementa- 
tion, we neglect star formation and therefore the energy feedback 
from supernovae and the gas mass drop-out term related to the star- 
forming phase is missing, leading to a higher baryon fraction of 
the hot gas phased compared to more self-consistent recipes (e.g. 
iBorgani et al.l2008l and references therein). In our runs, the baryon 
fraction of the warm-hot medium is expected to be overestimated 
by a factor ~ 2 in the innermost cluster regions. Given that we 
excised the innermost 0.1i?2oo from the analysis of clumps of the 
previous section, this problem is not expected to be a major one. 

Our model that includes cosmic rays, is taken from 
IVazza et all j2012h . and com putes the run-time e ffects of diffu- 
sive shock acceleration (e.g. Kang & Jones 2007) in cosmologi- 
cal shocks: reduced thermalization efficiency at shock, injection of 
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Figure 9. Top panels: projected X-ray flux (fogioLx [erg/s]) of 4 resimulations of a ~ 3 • 1O 14 M0 employing different physical implementations (see main 
text for details). Bottom panels: projected spectroscopic-like temperature for the same runs (T in [K]). 



CRs energy characterized by the adiabatic index Y = 4/3, pressure 
feedback on the thermal gas. This model was developed to perform 
run-time modifications of the dynamical effect of accelerated rel- 
ativistic hadrons in large-scale structures. Given the acceleration 
efficiency of the assumed model dKang & Jonesll20070 the pressure 
ratio between CRs and gas is always small, ~ 0.05 — 0.1 inside 
R200 of our simulated clusters at z = 0. 

The feedback models discussed above are an idealized sub- 
set of feedback recipes that have been used by other authors (e.g . 
ICattaneo & Tevssier 2007; Sii acki et alJl2007l ; [Gaspari et al.l201ll) . 
For the goal of this paper, they just represent extreme cases of dra- 
matic cooling or of efficient heating in clusters, while the reality 
lies most likely in between. 

Figure [9] shows the projected X-ray maps and spectroscopic- 
like temperature maps for four re- simulations (a-b-d-f) for the 
merging cluster. 

The corresponding profiles of average clumping factor and gas 
density for this cluster (bottom panels) and for the relaxed cluster 
(top panels) are shown in Fig[l0l The trend with the physical im- 
plementations is quite clear in both cases. The inclusion of cooling 
causes a "cooling catastrophe" and the enhancement of the core gas 
density in both clusters. At the same time it also triggers the forma- 
tion of denser and brighter clumps at all radii. A weak amount of 
thermal feedback from a centrally-located AGN is found to reduce 
the gas clumping factor only by a moderate amount, and only very 
close to the cluster centre « 0.3i?2oo), whereas it has no effect 
at larger radii. In these two runs the luminosity of the brightest 
clumps can be 10 — 10 2 times larger than in non-radiative runs. On 
the other hand, with a more efficient feedback recipe (early diffuse 
pre-heating an moderate AGN thermal feedback at low redshift) it 
seems possible to quench the cooling catastrophe in both clusters 
centre, and to reduce at the same time the gas clumping factor at 



all radii, making the latter only slightly larger than in the simple 
non-radiative run. The inclusion of CR feedback at cosmological 
shocks is found not to affect our gas clumping statistics compared 
to non-radiative runs. This can be explained because in our model 
the strongest dynamical effect of CRs is found around accretion 
shocks, where the energy ratio between injected CR-energy and 
gas energy can be ~ 10 — 20 percent. However, filaments can en- 
ter far into the main cluster atmosphere avoiding strong accretion 
shocks, and thus keeping a smaller energy ratio of CRs. Given the 
strong role of filaments in enriching the ICM of clumpy material 
(as shown in Sec l3.ll ), the fact that the overall clumping statistics 
within the cluster is not affected by the injection of CRs at shocks 
is not surprising. 

Finally, we tested the effect of increasing spatial resolution on 
the gas clumping factor in our runs. On average, a larger resolu- 
tion in the non-radiative runs is found not to change the overall gas 
clumping factor within the cluster, although assessing its effect on 
the high density peaks of the clumps distribution is non trivial. In- 
deed, the increase in resolution may change the ram-pressure strip- 
ping history of accreted clumps, causing a s lightly different orbit 
of sub structures towards the end of runs fe.g. lRoediger & Briiggenl 
120071) . This may lead to time-dependent features in the gas clump- 
ing factor profile, i.e. the higher resolution run of the perturbed 
cluster shows a "spike" of enhanced clumping factor at ~ R200 
at z = 0, while no such feature is detected in the resimulation of 
the relaxed system. The large-scale average trend of the gas clump- 
ing factor in both systems, instead, presents no systematic trend 
with resolution, suggesting that overall the effect of resolution in 
non-radiative runs is not significant for the spatial range explored 
here. However, we caution that some of the parameters of bright 
gas clumps derived in this work should be subject to some evolu- 
tion with spatial resolution. The small number of clumps found in 
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Figure 10. Left panels: average profiles of clumping factor for a relaxed and a perturbed (i.e. with an ongoing merger) cluster in the mass range ~ 2 — 3 • 
1O 14 M0 resimulated with various physical recipes (plus one resimulation at higher resolution). Right panels: gas density profiles for the corresponding runs 
on the left panels. 



these two cluster runs is so small that we cannot perform meaning- 
ful clumps statistics as in the previous Section. However, we can 
already notice that the higher resolution increases the brightness of 
the clumps by a factor ~ 2 — 3. In addition, also the typical inner- 
most radius of bright clumps can be only poorly constrained by our 
runs, and based on Sec l3.2l our estimate of a typical < 50 kpc/h 
must be considered as an upper limits of their real size, because the 
FWHM of most of our clumps is close to our maximum resolution, 
and this parameter has likely not yet converged with resolution. 

In conclusion, the effect of radiative cooling has a dramatic 
impa ct on the proper ties of clumps in our simulations (as suggested 
bv lNagai&Laull201ll) . However, AGN feedback makes the X-ray 
flux from clumps very similar to the non-radiative case and the most 
important results of our analysis (Sec. 3. 1-3.2) should hold. How- 
ever, the trend with resolution is more difficult to estimate, given 



the large numerical cost of simulations at a much higher resolu- 
tion. We therefore must defer the study of gas clumping statistics at 
much higher resolution to the future. This will also help to under- 
stand the number of unresolved gas clumps that might contribute to 
the X-ray emission of nearby clusters. 



4 DISCUSSION AND CONCLUSIONS 

In this paper we used a set of high-resolution cosmological simula- 
tions of massive galaxy clusters to explore the statistics and effects 
of overdense gas substructures in the intra cluster medium. While 
the densest part of such gas substructures may be directly detected 
in X-ray against the smooth emission of the ICM, the moderate 
overdensity associated to substructure increases the gas clumping 
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factor of the ICM, and may produce a (not resolved) contribution 
to the X-ray emission profiles of galaxy clusters. 

We analysed the outputs of a sample of 20 massive systems 
(~ 10 15 M(T)) simulated at high resolut ion with the ENZO code 
dNorman et al.ll2007l : [Collins et alJl201Qh . For each object we ex- 
tracted the profile of the gas clumping factor within the full cluster 
volume, and within smaller volumes defined by projected sectors 
from the cluster centre. The gas clumping factor is found to in- 
crease with radius in all clusters, in agreement with iNagai & Laul 
fcOllh : it is consistent with 1 in the innermost cluster regions, and 
increases to C(/>) ~ 3 — 5 approaching the virial radius. Strongly 
perturbed systems (e.g. systems with an ongoing merger or post- 
merger systems) are on average characterized by a larger amount 
of gas clumping factor at all radii. This enhancement is associ- 
ated with massive accretion of gas/DM along filaments, which also 
produces large-scale asymmetries in the radial profiles of clusters. 
Obtaining an accurate estimate of the enclosed baryon fraction for 
systems with large-scale accretions can be difficult because of the 
significant bias of gas clumping factor in the derivation of gas mass, 
and in the departure of real profiles from the NFW-profile, which 
can bias the estimate of the underlying DM mass. In a realistic ob- 
servation of a narrow sector from the centre of a cluster, the esti- 
mate of the enclosed baryon fraction can be biased by ±10 percent 
in relaxed systems and by ±20 percent in systems with large-scale 
asymmetries (not necessarily associated with major mergers). 

In order to investigate the detectability of the high density 
peaks of the distribution of gas clumping factor, we produced maps 
of X-ray emission in the [0.5-2] keV energy band for our clusters 
at different epochs. We extracted the gas clumps present in the im- 
ages with a filtering technique, selecting all pixels brighter than 
twice the smoothed X-ray emission of the cluster. By compiling 
the luminosity functions and the distribution functions of the most 
important parameters of the clumps (such has their projected size, 
temperature and temperature contrast with the surrounding ICM) 
we studied the evolution of these distributions as a function of the 
mass and of the dynamical state of with host cluster, and of the 
epoch of observation. Summarizing our results on clumps statistics, 
we find that: a) there is a significant dependence on the number of 
bright clumps and the dynamical state of the host cluster, with the 
most perturbed systems hosting on average a factor ^2 — 10 more 
clumps at all radii and luminosities, compared to the most relaxed 
systems; b) the host cluster mass, on the other hand, does not af- 
fect the number of bright clumps (although the investigated mass 
range is not large, ^0.3 dex); c) about a half of detectable bright 
clumps is located in the radial range 0.6 ^ R/R200 ^ 1.2 from 
the projected cluster centre, while a very small amount of clumps 
per cluster (order 1) are located inside < O.6R/R200 ; d) the typical 
size of most of gas clumps (extrapolated by our data) is ^ 50 kpc/h. 
In a preliminary analysis based on a few re- simulations of two sim- 
ulated clusters with different recipes of CR-physics or AGN feed- 
back, we tested the stability of our results with respect to complex 
physical models. Radiative cooling dramatically increa ses the ob- 
servable amount of gas clumping ( Nagai & Laul 120 111) . However, 
the required energy release from AGN seems able to limit at the 
same time the cooling flow catastrophe, and the observed clumping 
factor down to the statistics of simpler non-radiative estimates. The 
increase of resolution does not significantly change the observed 
number of detectable clumps or the gas clumping factor of the ICM. 
However, resimulations at higher resolution produce a ~ 2 — 3 
larger maximum X-ray flux in the clumps, and therefore in order to 
assess the numerical convergence on this parameter further inves- 
tigations are likely to be necessary. Also the typical size of bright 



clumps cannot yet be robustly constrained by our data, which only 
allow us to place an upper limit of the order of ^ 50 kpc/h. 

We conclude by noting that, based on our results of the statis- 
tics of bright clumps (Sec l3.2l ) it seems likely that a given num- 
ber of them could already be present in existing X-ray observa- 
tions by XMM and by ROSAT (and, to a much lesser extent given 
the expected lower statistics, by Suzaku). However, it is very diffi- 
cult to distinguish bright gas clumps from more point-like sources 
(like galaxies or AGN), given that their expected size is very small 
(< 50 kpc/h ), typically close to the effective PSF of instruments. 
In a real observation, a catalogue of point-like sources is generated 
during the data-reduction, and excised from the image. 

In order to perform a first check of the our quantitative re- 
sults, we examined the cat alogue of point- so urces generated in a 
real ROSAT observations of lEckert etail J2012h . In Fig[TT]we show 
the luminosity and the radial distribution of point-like sources from 
the ROSAT observation of A2142, and the corresponding simu- 
lated statistics obtained by analysing the whole simulated dataset 
at the same resolution and flux threshold of the real observa- 
tion. The luminosity distributions are taken for two radial bins: 
0.1 < R/R200 < 0.6 and 0.6 < R/R200 < 1.2. 

The number of point-like sources in the real observation is 
~ 1.5 — 2 times larger than the simulated distribution of clumps at 
all radii. The observed luminosity distribution inside the innermost 
radial bin also shows an excess with respect to the simulated clumps 
at all luminosities. In the outer bin, however, the observed luminos- 
ity distribution has a totally different shape with respect to the simu- 
lated one. It should be noticed than in the outer regions the sensitiv- 
ity to point-like sources in the real observation is much degraded. 
Therefore, it appears likely that the change in the observed lumi- 
nosity distribution function of point-like sources is mainly driven 
by instrumental effects. We found very similar results in the cat- 
alog ue of point- sourc es of the ROSAT observation of PK0754-191 
ofEc kertetaD ( l2012h . Based on such small statistics, and given the 
degrading sensitivity of ROSAT to point-like sources at large radii, 
it is still premature to derive stronger conclusions. The number of 
expected point-like sources brighter than ~ 3 • 10 _14 erg/ (s • cm 2 ) 
in the field is ~ 20/ deg 2 , based on the collection of wide field and 
deep pencil surveys performed with ROSAT, Chandra and XMM by 
iMoretti et all (120031). Given th e projected area of the A2142 ob- 
servation of lEckert et al.l d2012h . the expected number of point-like 
field sources is ~ 6 — 8 inside R200, significantly smaller than 
both the point-like sources detected in the field of A2142 and the 
total number of simulated clumps (that are in the order of ~ 40 
within R200). However, one should note that the density of galax- 
ies in a cluster field is much highe r than in a typical field, and so 
is the number of AGNs (see e.g. iMartini et al. 200^; lHaines et al. 
120121) . Thus, the larger number of detected point sources compared 
to the expectation from background AGN cannot be directly used 
as evidence for gas clumping. 

It is possible that a significant fraction of such point-like 
sources are compact and bright self-gravitating gas clumps in clus- 
ters. However, it is unlikely that they can be identified based on 
their X-ray morphology, given their small size. In the case of high- 
resolution Chandra images of nearby clusters, it is possible that 
some existing data can actually contain indications of gas clumps. 
Their detection, however, is made complex by the strong luminos- 
ity contrast required for them to be detected against the bright in- 
nermost cluster atmosphere. Very recently, the application of so- 
phisticated analysis techniques has indeed suggested that some gas 
clumps (with size ^ 100 — 2 00 kpc) may actually be present in 
nearby Chandra observations dAndrade-Santos et"al]|2012l) . 
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In the next future, cross-correlating with the additional infor- 
mation of temperature in case of available X-ray spectra may unveil 
the presence of gas clumps in the X-ray images (Fig.7). 
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